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Abstract 

Quantum fluctuations of some systems vanish not only in the limit 
h — > 0, but also as some other parameters (such as jj, the inverse of 
the number of 'colors' of a Yang-Mills theory) vanish. These lead to 
new classical limits that are often much better approximations to the 
quantum theory. We describe two examples: the familiar Hartree- 
Fock-Thomas-Fermi methods of atomic physics as well as the limit of 
large spatial dimension. Then we present an approach of the Hecke 
operators on modular forms inspired by these ideas of quantum me- 
chanics. It explains in a simple way why the spectra of these operators 
tend to the spectrum of random matrices for large weight for the mod- 
ular forms. 
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1 Introduction 



It is well-known that a classical mechanical system has many possible 
quantizations. The classical theory is the limiting case as ft —* , so 
it is not surprising that there would be many quantum theories that 
in this limit reduce to the same classical theory. In this, largely ex- 
pository, paper I will discuss the opposite phenomenon: how the same 
quantum theory can be obtained by quantizing radically different clas- 
sical systems. Viewed another way, a quantum theory could depend 
on two parameters, say h, A and the quantum fluctuations of some 
class of observables are of order fi\.. Then both the limits H — > and 
A — > are classical theories. These classical theories could be entirely 
different. In an example from atomic physics, the conventional classi- 
cal limit has a finite number of degrees of freedom, while the new one 
has an infinite number. We will refer to the the new limits (obtained 
by taking parameters other than Ti to zero) as 'neo-classical limits'. 

This phenomenon is physically interesting because one of the new 
classical limits may be a better approximation to the quantum theory 
than the naive classical limit. These ideas came into the high energy 
physics literature from the work of 't Hooft and Witten on the large 
N limit of gauge theories. Witten Q in particular worked out several 
simpler cases to popularize the notion that even case N = 3 may be 
well approximated by the large N limit. But historically, the vari- 
ous mean field theories of condensed matter physics Q (the spherical 
model for example) and even the theory of Nuclear Magnetic Reso- 
nance can be thought of as precursors of these ideas. 

For example, in atomic physics, in the usual classical limit h ^ 
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there is no ground state: the hamiltonian is not bounded from be- 
low. However, the neo-classical limit (in this case a version of the 
Hartree-Fock approximation) has a ground state. Moreover, it gives 
an excellent first approximation to the ground state energy of the 
atom. The neo-semiclassical expansion gives a systematic way of cal- 
culating corrections to arbitrary accuracy, although the complexity 
grows rapidly with the desired accuracy. 

Another important example, discussed below, is also from atomic 
physics: the quantum fluctuations in electron distances become small 
in the limit as the dimension of space becomes large. We will derive a 
simple effective potential that explains the stability of the atom. For 
another approach to this see the work of Hershbach || . 

Some important corners of mathematics are also illuminated by 
this phenomenon. The theory of modular forms can be viewed as the 
quantization of a classical mechanical system whose phase space is the 
upper half plane. The limit as the weight of the modular form goes 
to infinity corresponds to a classical limit. But there is also another 
classical limit corresponding to letting the level (the area of the fun- 
damental domain ) go to infinity. These limits lead to interesting new 
approaches to the problem of determining the spectrum of the Hecke 
operators on modular forms. 

The mathematical formulation of a classical dynamical system has 
expanded steadily in generality throughout history: as new physical 
theories are discovered we are led to enlarge the formalism to incor- 
porate the new developments. In the progression from ordinary differ- 
ential equations to Hamilton-Jacobi theory, symplectic geometry and 
the currently fashionable Poisson algebra formulation, we learn to deal 
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with increasingly sophisticated systems and symmetries. Neo-classical 
limits produce classical systems of even greater generality, often with 
non-local action principles and no simple hamiltonian description 
We have described some examples of this before. There are classical 
limits of quantum field theories that retain asymptotically freedom 
and require a renormalization ||. The challenge of finding the right 
mathematical description of these new kinds of classical systems re- 
mains. 

Much of the story told in this paper is, of course, well-known. I 
hope that organizing them in this way will help to understand common 
themes in apparently distant subjects. 



2 Hartree— Fock Theory of Atoms 
2.1 The Classical Limit of the Atom 

We start with a basic problem of quantum mechanics that cannot be 
solved exactly: an atom (or ion) with more than one electron. We 
usually start with the classical hamiltonian 

m p 2 m Ze 2 e 2 
^ On 5^ \ r .\ \ r . ~.| (^") 

i=l ^ i=l I 7 * I l<i<j<m \' 1 T l\ 

As usual, Z is the atomic number of the nucleus, e the charge of the 
electron and \i its mass. For each 1 < i < m, the position r$ and 
momentum pi of the electron are vectors in the Euclidean space 
Then we pass to the quantum theory whose hamiltonian is 



2 m V| m Ze 2 ^2 ^ ( ) 

i=i 2 V i=i N i<i<j<m \ n ~ r i\ 
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This is an operator on the complex Hilbert space of anti-symmetric 
wavefunctions T m = A m (7i). The anti-symmetry incorporates the 
Pauli exclusion principle. The space of single particle wavefunctions 
is 7i = L 2 (i? 3 , C 2 ); the wavefunction of each electron takes values in 
C 2 , since it can exist in two spin states 0. 

The limit fi — > is the usual classical limit of the atom. It is well- 
known that this is a spectacularly bad approximation to the quantum 
theory of the atom. The quantum hamiltonian is self-adjoint and 
bounded below and hence has a well-defined ground state. Indeed 
the central problem of atomic physics is the determination of this 
ground state wavefunction and the corresponding eigenvalue. The 
hamiltonian of the classical limit on the other hand, has no ground 
state: we can let position of an electron approach the nucleus, r a — > 0, 
thus descreasing the energy down to — oo. 

What is missing here is the uncertainty principle: in the quantum 
theory it is not possible to make the position of the electron close 
to the nucleus without making its kinetic energy large. This might 
suggest that there is no way to produce a classical approximation to 
the atom with a stable ground state. 

We will now produce a completely different classical system (with 
an infinite number of degrees of freedom in fact) whose quantization 
yields exactly the above quantum theory of the atom. Moreover, it 
has a ground state which is even a good approximation to the quan- 
tum ground state. What we will describe is just a reformulation of the 
standard Hartree-Fock approximation in atomic physics. This refor- 
3 We will, for simplicity, ignore relativistic and spin-dependent terms in the hamilto- 
nian. 
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mulation allows a generalization to relativistic many-fermion systems 
which we have described elsewhere 0. 

2.2 The Neo-Classical Theory of the Atom 

There are many standard texts that discuss the material in this sec- 
tion, although often without the geometric interpretation in terms of 
the Grassmannian. See for example [Q] Let Ti = L 2 (R 3 ,C 2 ) be the 
familiar complex Hilbert space. We define the Grassmannian Gr m (Ti) 
to be the space of linear self-adjoint trace-class^] projection operators 
of rank m: 

Gr m (H) = {p : Ti -> Ti\p* = p; p 2 = p; tr p = m}. (3) 

It is clear that an eigenvalue of p is equal either to zero or to one. Cor- 
responding to each such projection operator there is a subspace of Ti 
of dimension m: the eigenspace of p with eigenvalue one. Conversely, 
each such m-dimensional subspace V defines an orthogonal decom- 
position Ti = V © V^; then we can construct p as the hermitean 
projection operator to V. Thus we can see that Gr m (Ti) is really the 
set of all m-dimensional subspaces of Ti. Thus we have an infinite 
dimensional (but finite rank) generalization of the usual definition of 
the Grassmannian H. Gr m (Ti) is an infinite dimensional manifold, 
whose tangent space is the space of rank m self-adjoint operators on 
Ti. 

A classical dynamical system is specified by (i) a manifold which 
will be its phase space, (ii) a symplectic form on this phase space 
If a projection operator is trace class, its trace must be an integer, the dimension of 
the vector space it projects. 
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which will determine the Poisson brackets, and (hi) a real function on 
the phase space which is its hamiltonian. 

In our theory, Gr m (TC) is the phase space. The symplectic form 
on it generalizes the standard symplectic form on finite-dimensional 
Grassmannians. The Poisson brackets that it implies for a pair of 
functions is: 

{/,<?}= tr p[df,dg]. (4) 

Here df is the infinitesimal variation of / which can be thought of as 
a linear operator on 7i. If we represent the operator p by its integral 
kernel p b (x,y) ( where a = 1, 2 labels spin) we can write this Poisson 
bracket as 

{p%(x, y),p c d (z, u)} = 5 c b 5(y, z)p a d (x, u) - 5 a d 5{ Xl u)p c b (z, y). (5) 



The last piece of information is the hamiltonian, which we postu- 
late to be 

Ze 2 



#1 = /^^(^3-/^(^^ 3 
1 



+ - 



\x - y\ 



p a a (x,x)p b b (y,y) - p a b (x,y)p b a (y,x) d 3 xd 3 y{6) 



Here, p is the symbol of the operator p: 

pt(x,y) = 



pKx+^x-^e-^d'u, 



pt(^V i ,P)e^ ix - y) 



d 3 p 



(7) 



Clearly, p a a {x,x) = f ft(x,p)j0p. 

So far it is clear that this system depends on the parameters 
p, e, Z, m of the of the atom. Although the dynamical variables are 



operators, and h appears in the formula for the hamiltonian, it is a 
bona fide classical dynamical system. 

We will show that a quantization of this system is exactly the quan- 
tum theory of the atom. The physical meaning of the operator p is that 
it is the 'density matrix' of the electrons. Indeed p%(x, x) is the num- 
ber density of the electrons at the point x; tr p = J p®(x,x)d?x = m 
is the total number of electrons. More generally, p®(x,p) is the density 
of electrons of momentum p and position x. The Pauli exclusion prin- 
ciple which allows for at most one electron per single particle state, 
becomes the condition that this density matrix be a projection opera- 
tor, so that its eigenvalues can only be zero or one. Thus this classical 
dynamical system realizes many of the facts we usually associate with 
quantum theory. 

The first term represents the kinetic energy and the second term 
the potential energy due to the nucleus. We can combine these 'single- 
particle' terms in the hamiltonian into the form 

trpK, K = -^-V 2 -^. (8) 



Now, K is bounded below by E\ = — \p z ^ , as we know from 
the elementary theory of an ion with one electron. Since p is positive 
and tr p = m, we see that tr pK > mE\. ( A stricter bound can 
be obtained using the fact that p is a projection. But we don't need 
it.) This way, the system avoids the catastrophe of the conventional 
classical limit. 0. 

5 Strictly speaking H\ exists only on some dense domain of Gr m (H). This domain 
should be some class of pseudo-differential operators, and is the true phase space of the 
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The interaction of the electrons induces two kinds of terms. The 
first is obvious, the Coulomb energy of a charge cloud of density 
ep®(x,x). The last term is not so obvious-it is the 'exchange energy'. 
It is needed to get back the correct quantum theory (see below). By 
a version of the Schwarz inequality it should be possible to see that 



/ 



e 
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p a a (x, x)p b b (y, y) - p a b (x, y)p b a {y, x) ^--d 3 xd 3 y > 0. (9) 

x - y\ 

This expresses the physical fact that the electron-electron interaction 
is repulsive. Thus the total hamiltonian is bounded below by at least 
mE\. 

To actually find the ground state of this classical system, we must 
vary the hamiltonian subject to the constraints on p. Such a variation 
of p is always of the form 5p = —i[p,u] for some hermitean operator 
u. The condition for an extremum is then 

[p,dH 1 } = 0. (10) 

Here dH\ = is a linear operator 

dH 1 = K + U + W. (11) 

Here, K is as defined above and U is the multiplication by the 'mean 
field' 

r Zp 2 

u b( x ) = S b] J^PaiViV^y- (12) 
The 'exchange energy' contributes an operator W whose integral ker- 
nel is 

WZ(x,y) = - T ^-pt(x,y). (13) 
\x-y\ 



system. The correct statement is that the hamiltonian is bounded below within this 
domain. This is a technical project that I am unable to complete. It would be interesting 
to produce a functional analytic realization of the physical ideas that are described here. 
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Thus we can find the extremum by solving a non-linear eigenvalue 
problem self-consistently. We have, 

m 

P = X>a ® ^ (14) 
a=l 

where each of the vectors if) a £ 7i is an eigenstate of dH\. 

This is exactly the Hartree-Fock approximation to the atomic ground 
state. We find the wave-functions that are eigenstates of some sin- 
gle particle hamiltonian; the potential in this hamiltonian is self- 
consistently determined by postulating that m of these are occupied 
by electrons. Our description avoids the usual Slater determinants for 
the wavefunction-the hamiltonian only depends on the density matrix 
of the electrons and not the wavefunction itself. We have shown that 
this way of formulating the Hartree-Fock theory allows for general- 
ization to systems containing an infinite number of fermions such as 
relativistic theories [||. ||. 

Our point in this paper is that this theory can be thought of as 
minimizing the hamiltonian of a classical system on Gr TO (TC). This 
extends to the time evolution as well: the hamiltons equations of our 
system are the usual equations of time-dependent Hartree-Fock theory. 

2.3 Back To the Quantum Theory 

How do we quantize a system whose phase space is a Grassmannian? 
It is not possible to cover the Grassmannian by a single co-ordinate 
system, so it is inconvenient to look for canonical variables. However, 
the Grassmannian is a Kahler manifold, and we can apply the ideas of 
geometric (or Berezin-Toeplitz) quantization ||. In an earlier paper 
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(the appendix of ||) we used the representation theory of the unitary 
group to quantize this theory. 

Recall the situation in the case of finite dimensional Grassmanni- 
ans: let V be a finite dimensional vector space, and Gr m (V) the set of 
its m-dimensional subspaces. Gr m (V) is a compact Kahler manifold. 
Its canonical line bundle C admits a hermitean metric and a connec- 
tion whose curvature is just the symplectic form. ( A line bundle that 
admits such a metric and connection is said to be quantizable ]|] .) The 
holomorphic sections Hol(£) of this line bundle form a finite dimen- 
sional vector space isomorphic to A m (V). This space of holomorphic 
sections is a subspace of the Hilbert space of square-integrable sections 
with a projection operator II : L 2 (C) — > Hol(£). These geometric facts 
can be used to construct a quantization || of the dynamical system 
whose phase space is Gr m (V). 

From any function / : Gr m (V) — > R we will construct an operator 
/ : Hol(£) -> Hol(£) by the formula 

/ = n/. (15) 

That is, we multiply a holomorphic section by the function to get 
a section of C that may not be holomorphic; then we simply project 
out the holomorphic part. The operator we construct this way is self- 
adjoint (it is just a finite dimensional hermitean matrix in fact). 

In what sense is it a quantization of the dynamical system on 
Gr m (V)? How will we recover the classical limit? The idea || is that it 
is merely a special case of a one-parameter family of quantum theories 
where C above is replaced by C . As long as N is a positive integer 
the above ideas go through: there is still a projection Ily : L 2 (C N ) — > 
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Hol(£ N ) to a finite dimensional space of holomorphic sections. [] Also 
in the limit N — > oo the operator algebra tends to the Poisson algebra 
of functions in the sense that 

\\n N fii N g-ii N (fg)\\ = o(±), (16) 

and moreover 

\\iN[U N f, U N g]-U N ({f,g})\\=0(^). (17) 

Also the operator norm of T/v/ approaches the sup norm of the func- 
tion /. 

These ideas also extend to Gr m (7i) when V is replaced by the 
infinite dimensional vector space H. The technical aspects are simpler 
than in pj]l , since we need only finite rank projections. We give 
only a very brief outline here. Any subspace of dimension m can 
be brought to some standard subspace whose orthogonal complement 
is 7i±; hence Gr m (Ti) is a coset spaceQ Gr m (Ti) = Uo(H)/U(Tl±) x 
U(m). Using the trivial representation of U(7i±) and the determinant 
representation of U(m), we can construct a line bundle 

C = (U (H) x C) /U(H±) x U{m). (18) 

The holomorphic sections of this bundle can now be constructed and 

shown to form A m (7i). It is thus clear that A m (7i.) is the Hilbert space 

6 The space Ho\(£ N ) carries a representation of the unitary group Z7(V) given by 

the Young diagram of height n and width N, generalizing the completely anti-symmetric 

tensor representation of the case N = 1 above. 

7 I/o(7i) is the group of unitary transformations that mixes the standard m-dimensional 

subspace with its orthogonal complement only by a finite rank operator. It is analogous 

to the 'restricted Grassmannian' of Sato j^J except that it is modelled on finite rank 

operators rather than compact ones. 
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of at least one way of quantizing our system on Gr m (H). Indeed, the 
hamiltonian of the system when worked out in this way is exactly the 
quantum hamiltonian of the atom we had earlier. 

This quantum hamiltonian is the special case as N = 1 of a one- 
parameter family of theories. For N > 1 these describe fermions that 
carry a 'color' quantum number, except that only observables that are 
invariant under U(N) are realized in the Hilbert space Hol(£ N ). The 
Hartree-Fock method thus approximates the theory for N = 1 by the 
neoclassical limit as N — > oo. In effect -k measures the size of the 
quantum corrections. 

A relativistic generalization of this theory is described in § . There 
I developed an approach to two-dimensional QCD where the large N 
limit was realized as a classical theory. The story above appears in 
the appendix to that paper. Later I found that many other problems 
in physics and mathematics can be thought of in a unifying way as 
different classical limits of the same quantum theory. 



3 The Thomas— Fermi Approximation 

The problem of minimizing the energy on the Grassmannian is still a 
hard problem. Further approximations are needed. It turns out that 
there is a way to consider the limit[] h — » (a kind of semi-classical 
approximation) which yields a simpler theory. The ideas go back to 
the Thomas- Fermi approximation of early atomic physics Jill and have 



8 It is important to take Ti — > 0, after the theory has been formulated on the Grass- 
mannian as above; if we took h — > we would get a theory without a ground state, as 
explained above. 
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seen several revivals. There seems to be a connection with the density 
functional [14] method as well. Our point of view is based on symbol 



calculus and was in part inspired by the work of Lieb, Thirring [15] 
and others on the stability of matter. We will work out explicitly the 
leading terms but indicate how higher order terms can be calculated 
systematically if needed. 

It is possible to develop the theory in a more general context than 
in the last section without much additional workP| We therefore con- 
sider a system of m fermions with the hamiltonian 

m 

H = Et T (-^) + U ( r i)] + E G(Ti, rj ). (19) 

i=l i<j 

The configuration space of each fermion is R n . We will allow the 
fermions to carry a 'spin' quantum number a = 1, • • • Nf. The above 
hamiltonian is assumed to be independent of this quantum number. 
(We will usually suppress the spin index). Here, T(p) is the 'dispersion 
relation'; i.e., the dependence of kinetic energy on momentum. T(p) is 
usually spherically symmetric. U(x) is the external potential that all 
the fermions are subject to and G(x, y) = G(y, x) is the two body 
potential. 

In the last section we had the following special case: the dimension 
of space n = 3, the spin takes takes two values, Nf = 2, the kinetic 

2 

energy is T(p) = ^- and the inter-electron potential is the Coulomb 

2 

potential G(x,y) = \%_ y \ ■ The cases n = 1,2,3 for the dimension of 
space are also of interest in other contexts. 



9 In fact I worked this theory in this way (in 1992) to find a relativistic generalization 
of the Thomas-Fermi method. It remains unpublished. 
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We have then the Hartree-Fock energy 

H x (p) = Tt(T + U) P + 

-J dxdyG(x,y)[ tr p(x,x) tr p(y,y) - tr p(x, y)p(y,(@ff) 

to be minimized over all operators satisfying p^ = p, p 2 = p, Trp = m. 
We denote by tr the trace over flavor while Tr includes the integral 
over position as well. As before, the first term represents the single- 
particle kinetic energy and potential energy, the second term the direct 
interaction and the last term the exchange interaction. 

The minimization problem above leads to the variational equations 

[p,dffi]=0 (21) 

where the Hartree-Fock self consistent hamiltonian itself depends on 
p. A self-consistent solution is clearlyf^ p = Q(Ep — Hi) where the 
'Fermi energy' Ep is determined by the condition Trp = m. It is 
often too hard to solve this problem, so yet another approximation 
is needed. We could minimize over some smaller set of operators p 
thereby obtaining a variational bound that is simpler to calculate. Or 
we could calculate the function H\(p) in a semi-classical approxima- 
tion. 

The essence of the Thomas-Fermi approximation (a modern ver- 
sion is the density functional method pip is a combination of these 
two ideas: 

1. use the variational ansatz p = 0(— h) where h = t(—id) + v{x) is a 
separable hamiltonian (i.e., a function of p alone plus a function of x 
10 The theta- function of a self-adjoint operator Q(A) is denned to be the projection 
operator to the subspacc on which A is positive. 



16 



alone) ; 

2. expand the energy function 



H^ei-t-v)) 



(22) 



semiclassically 

3. minimize the leading term in this expansion HrFit, v) with respect 
to the variational parameters t and v. ( In many treatments, however, 
t is chosen to be the same as T(p) and only v(x) is varied.) 

The semiclassical expansion will amount to an expansion in pow- 
ers of the derivatives of v. The above ansatz for h is motivated by 
the form of the Hartree-Fock hamiltonian. If the two body potential 
G is absent, the first step is automatic, since H\ is already in this 
form. Even for interacting fermions, the direct interaction is already 
of the separable form. The indirect energy may not be separable in 
general, but as long as it is a monotonic function of momentum, the 
projection operator Q(Ep — Hi) will agree with that of some separa- 
ble hamiltonian. Thus one expects this separable ansatz to be a good 
approximation. 

The projection operator Q(h — E) can be expressed in terms the 
resolvent operator, i-zt?, since 



where D is a contour that surrounds the negative real axis in a coun- 
terclockwise direction. There is a semi-classical expansion for the re- 
solvent which can be used to derive one for the projection operator 




(23) 



e(h-E). 
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3.1 The Semiclassical Expansion of the Re- 
solvent 

To do the semiclassical expansion, it is convenient to restate the prob- 
lem in terms of (matrix-valued) functions on the phase space R n © R n 
rather than operators on the Hilbert space L 2 (R n ,C N f). There is a 
systematic theory of this procedure (symbol calculus) decribed in de- 
tail in for example, Ref. The main idea is to use Weyl ordering 
to set up a one-one correspondence between functions on the phase 
space and operators on the Hilbert space. From a function A(x,p) we 
construct the operator A whose kernel is 

A(x,y) = J A^±y,p)e L n^-y)[dp\. (24) 

We use the abbreviation [dp] = ^ 2 th) n • ^ we a PPh/ this to simple 
functions such as polynomials we can check that this definition corre- 
sponds to Weyl ordering. For example, the function xp becomes the 
operator —ihdx + x(—ihd). 

Conversely, given an operator, we define its symbol to be the func- 
tion 

A(x,p) = J A{x + ^x- ^)e~y- z dz. (25) 

The idea is that p is the momentum conjugate to the relative coor- 
dinate of the operator kernel. The operator multiplication can now 
be translated into the multiplication of symbols. The result can be 
expressed in closed form: 

AoB(x,p) = L~ i ^^~ WIIJr) A(x,p)B(x',p')\ . (26) 

^ J x=x';p=p' 



18 



The trace of operators becomes an integral in phase space 

TrA = tr J A(x,x)dx= tr J dx[dp]A(x,p). (27) 

We emphasize that the algebra of symbols under the multiplication 
law is exactly the same (isomorphic) to the algebra of operators on a 
Hilbert space;i.e., no approximation is involved in replacing an oper- 
ator by its symbol. 

We see that to the leading order the above muliplication law is 
just the pointwise multiplication of the classical theory In the next 
order there is a correction proportional to the Poisson bracket. If we 
expand the exponential, 

AoB(x,p) = Y / { — ) -M,B} [n) . (28) 

?1=0 ^ ' 

Here we see a sequence of generalized Poisson brackets 

n 

iAM^ = Yy^^-Z-M^r (29) 

r=0 

where A 1 = and A\ = etc. n = 1 corresponds to the usual 
Poisson bracket. If A and B commute as matrices on spin, the odd 
brackets are antisymmetric and the even ones are symmetric. Other- 
wise, there is no particular symmetry property. 

Consider now the resolvent operator of a hamiltonian h, 

r(E) = (30) 

We will now derive a semiclassical expansion for the symbol f(E) of 
this operator. The resolvent symbol satisfies 

r(E) o(h-E) = l. (31) 
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Expand r(E) in power series in % and put into the expansion of the 
above equation to get, 

oo 

r(E) = J2r {k )(E)h\ (32) 

00 00 / -jh\ n 1 

EE — ^*{f (Jfc) (^,fc- J B} (B) = l. (33) 

n=0 fe=0 

Equating the powers of h on both sides of this equation, we get a set 
of recursion relations 

f {0) (E) = (h-E)-\ (34) 

m ■ -1 

f {m) {E) = -Y J {^) n ^{r { n-rn){E)M(n){h-E)-\ (35) 

n=l 

If /i is diagonal in flavor space, the odd terms r/ 2m +i) vanish. The 
above expansion can be used to derive the usual WKB quantization 
conditions as well as higher order corrections to it. 

Of particular interest to us is the case where h is a separable op- 
erator: 

~h(x,p)=t(p)+v(x). (36) 

In this case the mixed derivatives in the generalized Poisson brackets 
vanish and we get 

{f {k) M(n) = r {k)h ... in t h - in + (-l) n ^p^... in . (37) 
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Iff"] t(p) = piPi as in nonrelativistic quantum mechanics, there is 
only one term for n > 2, 

{f(k)M(n) = {-l) n r^ ln v h ... ln (38) 

while 

{r(k),h}(l) = 2 Pif(k)i ~ r\ k) Vi (39) 

and 

{f(k)M(2) = 2 ^(Jfc)M + f (fe)%- ( 40 ) 

If moreover, v is diagonal in flavor space, fm = and 

f « (E) = i(^ l S + 2fe '"' 2 -"«>- (41) 



3.2 Derivative Expansion of Energy Function 

Now we can rewrite the Hartree-Fock energy in terms of the symbol 
p(x,p) of the projection operator p. 

H\(p) = tr J p(x,p)T(p)dx[dp] + tr J U(x)p(x)dx+ 
- J G(x,y) tr p(x) tr p(y)dxdy 
-\J dx J[dp][dp']G{x,p - p') tr p(x,p)p(x,p') (42) 

Here, 

p{x) = p(x,x) = J [dp]p(x,p). (43) 
We use units here such that 2\i = 1, to simplify the formulas. 
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We must minimize this subject to the constraints 

p o p(x,p) = p(x,p); tr J p(x,p)dx[dp] = m. (44) 

We reiterate that although the problem has been formulated on the 
classical phase space, no approximation has been made yet. All the 
complications are in the multiplication law of the functions (hence in 
the quadratic constraint on p). 

Now we put in the separable ansatz (which satisfies the constraint 
automatically) and expand in powers of H. We will have 

oo 

p{x,P) = P(k)(x,p)h k (45) 

fc=0 

Using the integral representation in terms of the resolvent symbol, 

dE 



/ —f(E) 
Jd 2m y ' 



(46) 
we get 

<W D 5§^> (£) ' (47) 

These terms in the expansion of p are distributions on the phase 
space involving the delta function and its derivatives, although the 
f(fc) are ordinary functions. In the same way, we have expansions for 
the number density, 

oo 

P {x) = Y,~P{k){x)fi\ (48) 

k=0 

2~ J[dp]f {k) (E,x,p) (49) 
and the kinetic energy K = tr / T(p)p(x,p)dx[dp], 

oo 

K = J2K {k )h k . (50) 

k=0 
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Also, 

f dE f 

Q(k) = J d ^— tr J dx[dp]r (k) (E,x,p)T(p). (51) 

The direct energy can be written in terms of the density function 
p(x). The exchange integral is more complicated, being quadratic in 
p; however, in most cases it is quite small and explicit calculation in 
higher orders is not necessary. 

It is now straightforward to calculate the Thomas-Fermi energy to 
lowest order in the case of nonrelativistic quantum mechanics with a 
potential v that is diagonal in flavor space. We get upon evaluating 
the integrals, (it is convenient to introduce a new variable by v(x) = 



n 



Here 



< = Wr-^W) (53) 

is the area of a sphere of unit radius in momentum space. The 
exchange integral is,to lowest order, 



1(0) 



= \l [dpdp']Q{<j ) \x)-p 2 )Q{ ( p 2 {x)-p' 2 )G(p-p'). (54) 



=2 



With G{p) = 3 as for the Coulomb interaction, we can evaluate this 
more explicitly by introducing spherical polar coordinates in momen- 
tum space. We get 



ho) = \^^C n I ^~\x)dx (55) 
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where 



i rl . ^ fir s in n - 2 > 



JO Jo 



C„ = / dy / tffaftW rfg 2 ^ /2 9 ; (56) 

y + y — zyy cos ft 



This leads to 
EtfW = tr/K^(*)-^^|^ 



\^\ G(x, y) tr 4>{x) n tr <t> n (y)dxdy. (57) 
2 n z J 

In our expansion the exchange term appears in the lowest order. How- 
ever, in atomic physics it is as small as the terms involving derivatives 
of cj), so it is often ignored in the lowest order treatments. Also, in 
many discussions, the energy is expressed as a function of the den- 
sity p(x), but one can make the change of variable from the Fermi 
momentum 4>(x) to p(x) easily. 

Now we can vary this w.r.t. to eft to get an integral equation that 
determines the ground state in this approximation. Actually a more 
convenient variable to use is the mean field induced by this electron 
density: in terms of it we get a differential equation instead. If the 
distribution is spherically symmetric, as for an atom, this becomes a 
second order non-linear ordinary differential equation, the celebrated 



Thomas- Fermi differential equation [13]. 

Thus it is indeed possible to take the limit as h — > on systems 
such as the atom and get a sensible approximation to the ground state. 
However, this leads to a density function in the classical phase space 
and not conventional classical mechanics. Moreover, it has to be de- 
rived through an intermediary that is a bona-fide classical mechanical 
system but of infinite dimensions. 
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4 Atoms in the Limit of Large Dimen- 
sion 

As another example of a neoclassical limit, again in atomic physics, we 
consider the limit of large spatial dimension. This idea originates in 
an observation of Witten that in this limit the quantum fluctuations 
in the rotation invariant quantities will become small. Let r a i for 
a = 1, ■ ■ ■ m and i = 1, • • • n be the positions of m electrons in an atom 
(or ion) of atomic number Z. Although the physically interesting case 
is n = 3 we can, as a mathematical device, extend the system to 
n spatial dimensions. The problem of determining the ground state 
becomes that of minimizing 

n^-( 58 ) 

ai 

subject to the condition that 

J\iP(r)\ 2 Y[dr ai = l. (59) 

ai 

Here if; £ A m (L 2 (K n ,C N f)}. 

Now, the hamiltonian is invariant under the rotation group 0(n) x 
U(Nf). We now take the limit as n and Nf tend to infinity, and 
recover a classical theory. When Nf is large, we can assume that 
the wavefunction is completely anti-symmetric in the 'spin' indices; 
the position dependent part of the wavefunction is then symmetric. 
Indeed we can assume that this part is rotation invariant, so that 
it depends only the invariant quantities q ao = ^T^rfo. 
12 In the real world, the ground state wavefunction is invariant under O(n) at least for 
the noble gases: all the shells are filled. 
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h 2 di(>* di; / Ze 2 x ^ e 2 , . 
zjior ai or ai i a \r a \ 1<a<b<m \ r a r b\ 



These inner products form a positive m x m matrix. A com- 
plete set of 0(n) invariants are given by the remaining bilinears L b a = 

hV*3, im-\ + > pab = -fdS^i- They form a representation of the 
symplectic Lie algebra Sp(2n): 

[q ab ,qcd] =0= [P a \P cd ] 

[L%, L c d \ = — (s b L d — 5 d Lfj 

[L a b , P cd ] = ^ (5 c b P ad + 5 d P ac ) 

[L a b ,Qcd] -^(.S a c q bd + 5 a d q b c). (60) 

These commutators are proportional to -. Hence, there are two 
limits where the quantum fluctuations vanish: the conventional clas- 
sical limit where we let Ti — ► keeping n fixed (at the value 3 for 
example), or the neo-classical limit where we let n — > oo keeping Ti 
fixed. In this neo-classical limit, the quantum observables tend to 
classical ones satisfying the Poisson brackets of the sympleclassical 
theory.ctic Lie algebra: 

{qab,q cd } =0= {p- h ,P cd } 

{L a b ,L d } = h{5iL a d -5 a d Ll) 

{L a b , P cd } = h (5 c b P ad + 5tP ac ) 

{Llq cd } -h{5 a c q bd + 5 a d q bc ). (61) 

These Poisson brackets will determine the neo-classical equations of 
motion, once the hamiltonian is determined. 

There are some subtleties in determining the hamiltonian of this 
neo-classical theory: there is a new term in the potential arising from 
the change of the measure of integration. (It is possible to interpret 
this as a kind of 'Fischer information' while the measure determines 
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a kind of 'entropy' [|17[ . But we don't need this idea here.). Once the 
correct hamiltonian has been determined, this classical theory gives 
a relatively simple minimization problem for the ground state energy. 
Here we will consider only the static limit ( time independent solution) 
that determines the ground state of theory. 



4.1 The Change of Variables 

Let us return to the variational problem of determining the ground 
state. This will reduce to the minimization of an effective potential 
that depends only on q a t>. First of all, we need to determine the 
measure of integration fi(q) Ilc<d dqcd '■= n(q)dq determined by the 
change of variables q a b = ^r a iru on the Lebesgue measure Y\ a idr a i. 
This can be done by evaluating the following integral in two different 
ways: 

Z(J) = I e-^ jab fi(q) J] dq cd = f e^^" ]Jdr CJ . (62) 
J 9>o c <d J cj 

( Here, J is a positive matrix.) On the r.h.s. we have a standard 
Gaussian integral yielding 




e~ q « bJab ^q) J] d qcd = k(n, m) (det J)~^ (63) 



where k(n,m) = (Tin)~~ is independent of J. 

Thus Z(J) depends on J only through its determinant. It fol- 
lows0 that [i(q) can only depend on q through detg, prompting the 
13 The space of positive matrices is a homogenous space of the general linear group, 
since any such matrix can be mapped to the identity by the transformation q i— > SqS T 
with J transforming dually. The transformation law of Z{J) under this transformation 
completely determines that of fi{q) as well. 
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ansatz fi(q) = k[detq] u . To determine v we note that under the 
transformation q — ► SqS T , the measure of integrtion transforms as 
dq^ [det S] m+1 dq. Thus 

Z(J) = f e- trq ( sTjS 1 [det q}" [det S] 2u+m+1 dq = Z{S T J S)[det S] 2u+m+1 ,{U) 

Jq>0 

which determines v = n ~™ _1 . 
Thus we have 

M\ 2 = ~k[ \4>{q)[ 2 [detq] : ^ 1 dq. (65) 

Jq>0 

It is thus tempting to define a new wavefunction absorbing the de- 
terminant pf q: 

x(?) = V/*(#(?), IM| 2 =/ \x(q)\ 2 dq. (66) 

Jg>0 

This x{q) is a kind of 'radial wavefunction'. 



4.2 The Effective Potential 

Now we must express the hamiltonian in terms of this %■ The only- 
calculation we need is for the gradient of the wavefunction: 



/ 



h 2 dip* dip 
2/in 2 dr ai dr t 



Y[dr bj = 



h 2 i 

2/ a?- 6 



^x*)^djgx)^ 



2/i 5 
<9x 



ab cd 



+ 



Q q cd 

n — m — 1 d log det q 



dq 



cd. 



+ 



Q q ab ' 4 Q q 

n — m — 1 d log det q 



ab 



dq- 



cd 



x(q) 



dq (67) 



where 



ab cd 



n 



, dq ab dq cd 
dr ej dr ej 



S ac q bd + 5 aa q oc + goc^ + ^ 



sad be 



be ad 



Ad ac 



(68) 
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is an induced metric on the new configuration space. Moreover we 
know from eiementary matrix theory that 

dlogdetg _ x 

dq ab = lab ■ (69) 

The terms prportional to \*X become a correction to the potentiai; 
the terms invoiving one derivative of the wavefunction combine to give 
a total derivative that can be dropped. Those that involve the square 
of the derivative of \ become a new kinetic energy term. 
Thus the variational problem is now to minimize 



ft 2 ate**** d X 
2 fin 2 dq ab dq cd 



lq>0 

subject to the constraint 



^I^iSr + (Veff(q) + V(q)) \ X {q)\ 2 



dq (70) 



/ \x{q)\ 2 dq = l. (71) 

Jq>0 



Here, 



ft 2 (n - m - l) 2 ,. 

Also, U (q) is the potential energy of the electron expressed in terms 
of the new variables: 

u{q) = " S 7^ + v^ + J*-***] (73) 

where a = 4-. 

2 

We are now ready to take the limit as n — > oo, holding = a 
( not e 2 itself!) fixed. As expected, in the new variables, the kinetic 
energy of the ground state wavefunction will be of order It is 
very important that there is now a new term in the potential energy 
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( arising from the kinetic energy of the old picture) which makes it 
bounded below: 

V(q) = ^-tiq- 1 + U(q). (74) 

In the end the correction to the potential is quite simple!. 

The ground state energy in our neoclassical approximation is the 
minimum this function over all positive q. The condition for this is an 
algebraic equation for q. 

The case of a hydrogenic ion is of course simplest: when m = 1, q is 
just a positive number and there is no repulsive Coulomb interaction: 
-r j / \ h 2 _i Za 

V ® = 8? ~ ~q (75) 



The minimum is 



Z 2 a 2 2 Z 2 e 4 
2p-^- = --M-^-- (76) 



This is to be compared with the exact answer ( for n = 3) of — ^ z § . 
Thus we get roughly the correct answer: the relative error is about ^. 
It should be possible to improve on this by semi-classical methods. 

More generally, it is reasonable to expect (but not guaranteed) that 
the minimum will respect the permutation symmetry of the problem. 
Then we can put the ansatz that all the diagonal elements are equal 
(say, q aa = p 2 , Va) and that all the off-diagonal elements are also equal, 
(put q a i) = p 2 u,\/a b). Then |u| < 1 by Schwarz inequality. The 
potential becomes in these new variables, (it is convenient to choose a 
kind of atomic units during such explicit calculations, 2p = h = a = 

1): 



V(p,u) = ^f(u)-- p9 (u). (77) 
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Here, 



(78) 



Moreover, 

f(u)=trq~ 1 , q=(l-u)+uC (79) 

and C is the mx m matrix all of whose matrix elements are equal to 
one. The spectrum of C is quite simple: it has an eigenvalue equal to 
zero with degeneracy m — 1 and the remaining eigenvalue is just m. 
Thus we can determine the spectra of q and g _1 and hence its trace: 

771 — 1 i 

f(u) = ^ + - -. — . (80) 

JK ' l-u \ + {m-l)u y ' 



It is simple to minimize in p to reduce the problem to minimizing 



in u of 



9 2 (u) 



. If we change variables yet again to 



1 



1 

vW-u)]' 2-^' 



(81) 



our approximation of the ground state energy becomes the minimum 
of the rational function 

2mv 2 — (m — 1) 



V(v) 



m(m — 1) 1 2 
mZ v 



2mv 2 [2(m - l)v 2 - (m - 2)] ' 



(82) 



This minimum can in fact be found in closed form as an algebraic 
function of m and Z. But the formula (obtained by an algebraic 
computation program such as Mathematica) is quite complicated. But 
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this formula is fit very weliP| in the case of a neutral atom (i.e., m = Z) 
by the polynomial 

r i 2ue 4 

E(Z) = - 0.00152507 + 0.0871987Z - 0.920957Z + 1.83211Z 2 -^(83) 

We have restored the original units to make comparison with other 
methods easier. 

The point of this method is that it gives an exactly solvable and 
reasonably accurate picture for the ground state of the atom without 
having to deal with complicated nonlinear differential equations. The 
answers are reasonable considering the simplicity of the calculations. 



5 A Physicist's View of Modular Forms 

Next we will consider an example from mathematics: the theory of 
modular forms. I don't claim to have solved any deep problem in 
this area ( of which there are many). But perhaps the point of view 
described will suggest new methods. 



5.1 The Modular Group and its Subgroups 

We will give only a foretaste of the theory of modular forms. See Ref. 



[18] for most of the proofs and precise statements of the results. Also 



see Ref. [19| for relations to other areas of mathematics and physics. 

The group of two by two matrices with integer entries and de- 
terminant one is called SL2(Z); its quotient by the center, T(l) = 



14 



The fit is good to a relative error of 0.01% over the range 1 < Z < 100 
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SL2(Z)/Z2, is the modular group. It is conventional to denote ele- 

fa b\ 

ments of F as matrices , their pre-images in SL2{Z). 

\c <l ' 

r(l) acts on the upper half of the complex plane U through the 



fractional linear transformations 

az + b 



cz + d 

A fundamental region for this action is, 



(84) 



D = {z\\z\ > 1, |Re z\ < -}. (85) 

This region is a spherical triangle with vertices at ioo, ±i + ^ri. The 
point is that translations can be used to bring any point inside the 
strip | Re z\ < i; and under inversion any point is equivalent to one 
outside the unit circle. 

The modular group is generated by S : z i— > — ~ and P = ST : 
z i — > 7T. It is obvious that S 2 = 1,P 3 = 1. Indeed it can be 

z+l ' 

shown that T(l) = Z2 * Z3 is the free product generated by these two 
elements-there are no other relations among these generators. 

Many interesting groups appear as subgroups of the modular group. 
For example, the commutator subgroup of T(l) is the free group on 
two generators; it is a normal subgroup of indexF^J 6. Thus the mod- 
ular group is both non-abelian and infinite in an essential way: free 
groups are the ultimate examples of such groups. 

The principal congruence subgroup T(n) of level n consists of all 
elements that are equal to the identity matrix modulo n. Now we 



15 The index [G : H] of a subgroup H of a group G is the number of elements in the 
coset G/H; alternatively, it is the number of copies of the fundamental region of G that 
is needed to form a fundamental region of H . 
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see why the modular group is called T(l); its elements are of the 
(a b\ 

form with a = d = 1 mod n, c = b = Omod n. Any 

\c dj 

subgroup T in between, T{n) C V C T(l), is called a congruence 
subgroup of level n. The congruence subgroups are all of finite in- 
dex. It is possible to show by a counting argument [|l8| that [r(l) : 

r(n)] = n 3 IlpirJl ~~ P~ 3 ]- Of particular importance [^| is the subgroup 

(fa b\ ) 

ro(n) = < G r(l),c = mod n >. The index can be shown 

to be@ [r(l) : r (n)] = nl[ p]n [l + p" 1 ]. 



5.2 Modular Forms 

A entire modular form of integer weight k associated to a subgroup 
r C r(l) is a holomorphic function on the upper half plane (including 
the point at P] zoo) satisfying 

(cz + d)- fc /(^±^) =/(z), for ^ J er. (86) 

It is called a cusp form if /(ioo) = 0. 

If the weight is even, we can think of a modular form as a covari- 
ant tensor of order k/2 ('form') on U/Y: the condition above is the 
statement that f{z) [dz] 2 is invariant under Y. 

The homolomorphic sections of the canonical line bundle on U ( in 
this case the cotangent bundle) C are entire functions f{z) on U such 
that f(z)dz is invariant under Y. Thus the modular forms of weight 



16 It is the modular forms of weight two with respect to this subgroup that appear in 

the Shimura-Taniyama conjecture. 

17 A function / is holomorphic at ioo if it has a convergent Fourier expansion f(z) — 

E(T fne 27Tlnz . Moreover, /(ioo) = fo- 
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k are simply holomorphic sections of £2 . If A: is odd these correspond 
to some 'spinors' on U/T. 

An example of a modular formPI of weight 2k is the Eisenstein 
series 

E (87) 

(m,ra)^(0,0) v ; 

It does not vanish at ioo: G2fc(ioo) = 2((2k). 
The most famous cusp form is 

00 

AO) = (2ir) 12 e inz J[ [1 - e 2 ™*] 12 . (88) 
n=l 

It is of weight 12. It is nonzero everywhere except for a simple zero 
at ioo. It appears in Ramanujan's theory of partitions of numbers. If 
we expand the product of the twelveth root of A (which is called the 
Dedekind ^-function) we can see that it is a generating function for 
partitions. The partitions of large numbers is given by the asymptotic 
behavior as Im z —* 0; this is an essential singularity of the function so 
at first this looks hopeless. However, the modular invariance relates 
the value of A at z to its value at — -; thus the behavior at ioo (which 
is trivial to determine) gives the behavior as Im z — * 0. Hardy and 
Ramanujan turned this rough stone of an idea into an exquisite jewel 
(further polished by Rademacher) , deriving an asymptotic formula for 
partitions of large numbers. 

Any modular form is a periodic function hence can be expanded in 
a Fourier series. These Fourier coefficients are of great interest. An ex- 
ample is the Ramanujan r-function, which are the Fourier coefficients 



IN 



If we don't specify T, we will be speaking of the modular group itself. 
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of the modular form above, 

oo 

A(z) = Y / r(n)e 2ninz . 
1 

A deep conjecture of Ramanujan (proved eventually by Deligne fol- 
lowing ideas of Grothendieck) was that 

\ T (p)\ < 2p%. (90) 

In the theory of partitions, this inequality gives a bound on the error 
term to the Hardy-Ramanujan asymptotic formula for partitions of 
large numbers. These error terms seem to oscillate erratically yet a 
bound on their magnitude follows from the above inequality. 

These erratic oscilations are related to yet another interesting phe- 
nomenon: if we define 2p~z cos 9(p) = r(p), the angles 9{p) seem to 
be distributed randomly according to the circular ensemble of random 
matrix theory. Indeed spectra of random matrices appear in many 
places in the theory of modular forms and related Dirichlet series (see 
the recent books We will seek a clarification of this phenomenon 

using ideas from quantum mechanics in the theory of Hecke operators. 

Let M/% be the vector space of entire modular forms of weight k 
and 5fc that of cusp forms. It is clear that M^M\ C Mk+i. For k > 
12, multiplication by A gives a linear map M^_i2 — > S^- Moreover, 
dim5fc = dimM^ — 1 since there is just one condition on the Fourier 
coefficients of a cusp form: that the zeroth one vanishes. It is possible 
to reduce the determination of the dimension of to small values of 



k using these facts; see [18| for details. 

With our definition of weight, there are no entire modular forms 
of odd weight. For k = there is just one entire modular form, the 
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constant. There are none for k = 2. For k = 4,6, 8, 10 the only entire 
modular forms are multiples of the Eisenstein series. 
For k > 12 and even, 



if k + 2 mod 12 1 
dimS fc = { f ^ J . (91) 

1 if k = 2 mod 12 I 



Thus for large k, the dimension grows linearly with weight. A way of 
understanding this is that the elements of Sk are holomorphic sections 

k 

of the line bundle £2; as k grows this line bundle has greater Chern 
character allowing for more sections: it approaches a kind of classical 
limit. 

5.3 Modular Forms as Wavefunctions 

To a physicist, the above theory of modular forms is very reminiscent 
of quantum mechanics. 

We can regard the upper half plane as the phase space of some clas- 
sical mechanical system. The symplectic form is the Poincare form: 

dx Ady 

lo = ^ — , z = x + iy. (92) 

y 2 

The modular group ( or one of its finite index subgroups) can be 
thought of a discrete gauge group, so that points related by such a 
transformation represent the same classical state. A wavefunction 
would be a holomorphic function on the upper half plane; more pre- 

k 

cisely it would be a holomorphic section of a line bundle £2 on U/T. 
Thus I is analogous to the parameter N in our earlier discussion of 
compact Kahler manifolds. The base U/T is not usually a compact 
manifold because of the cusps ( points at infinity and points where 
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the stability group is finite) . Nevertheless U /T has finite area hence 
most of the theory ought to generalize. 

What is the dynamical system whose phase space is U/T(l)7. We 
can imagine it as a model of quantum gravity in two dimensions. 
There are many such models that illustrate various aspects of grav- 
ity. Here, we think of space-time as a torus. Our model of gravity is 
conformally invariant (not crazy since two is the critical dimension for 
gravity). Thus the set of metrics modulo diffeomorphisms and con- 
formal ( Weyl) transformations is the phase space of gravity. Using 
a diffeomorphism that is connected to the identity and a Weyl trans- 
formation we can bring any metric to the form ds 2 = \d9\ + zdQi | 2 , 
where < 9i, 62 <2tt are standard co-ordinates on the torus. Also we 
can choose Imz > 0. Now if we also allow for diffeomorphisms that 
are not connected to the identity, which are 
f9 1 \ (a b\ /0A fa b\ 

then the true phase space would be U/T(l): the action of T(l) on z 
is exactly the above fractional linear transformation. 

What would be the meaning of a gauge group that is only a sub- 
group of T(l)? We might have some additional geometric object on 
the torus that has to be invariant as well ( like a spin structure ) that 
would reduce the size of the gauge group. 

What would be the hamiltonian of our theory? A closed cosmology 
like a torus would at first not seem to have any meaningful time evolu- 
tion. An asymptotically flat space-time would have a time at infinity 
with respect to which we can evolve its wavefunction. In closed uni- 
verses in four dimensions, the Wheeler-DeWitt equation gives a 'time 
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evolution' where the conformal factor of the metric itself is a kind of 
time variable. But we have given this up by postulating that contin- 
uous rescalings are part of the gauge group, so that the wavefunction 
is invariant under them. 

However we can still regard time evolution as a discrete rescaling ( 
'expansion') of the universe. While rescalings connected to the identity 
are part of the conformal group, there are certain discrete rescalings 
that for example double the size of the fundamental region. There 
are many different ways of rescaling such a region (e.g., double just 
one leg of the fundamental parallelogram) which individually violate 
modular invariance. Only by averaging over all of them would we 
recover modular invariance. 

What we describe above is an interpretation of the Hecke opera- 
tors on modular forms:they are rescalings averaged over the modular 
group. Because the evolution is discrete we cannot find a generator 
for infinitesimal transformations. The closest we get to are the prime 
rescalings, which cannot be decomposed as compositions of others. 
They yield a family of commuting hermitean operators which together 
play the role of the hamiltonian. In the next section we give a more 
detailed description of Hecke operators. 

The central problem of the Hecke theory of modular forms-determining 
the simultaneous eigenvectors of the Hecke operators- is just like the 
central problem in quantum mechanics- finding the eigenfunctions of 
the hamiltonian. Quantum mechanics suggests some strategies to at- 
tack this Hecke problem. The limits of large weight or large level are 
like neo-classical and classical limits. For example, the number of lin- 
early independent modular forms is O (uk) in the limit of either large 
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index v or large weight k. The number of independent states of a 
quantum mechanical system with a two-dimensional phase space is of 
order of the area of the phase space divided by h. The fundamental 
region D of the modular group is not compact but still has finite area 
with respect to the Poincare metric: 



The fundamental region of a subgroup of index v is then just uA(D). 
Since the number of linearly independent modular forms is | for large 
k, we see that the analogue of fi in our theory is essentially 6A(D)/k = 
t£ . For example for the subgroup Tq (p) which has index ~ p (for prime 
p) there should be, according to this interpretation, ~ pk linearly 
independent modular forms. This is indeed known in the traditional 
theory of modular forms. 

Thus we get simpler classical analogues of the theory of modular 
forms in these limits; we can then hope to understand the general case 
by asymptotic expansions in inverse powers of k or n. In the limit of 
large n, the subgroup 
fa b\ 

To(n) = { |c = Omod n; ad — be = 1; a, b,c,d G Z} (95) 

\c dj 

becomes essentially the group of translations 

( \;c£Z, z^z + c. (96) 

Vo \) 

This is a huge simplification: the invariance group becomes more 
'abelian' as n — > oo. In ordinary gauge theories ( such as Yang-Mills 
theories) the limit as the theory becomes abelian and the limit of small 
quantum corrections are intimately related. ( Perturbation theory is 
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essentially the same as the loop expansion.) Thus one elementary 
strategy to understand modular forms is to study first this easy lim- 
iting case where modular forms reduce to periodic functions on the 
upper half plane. We will see then that this 'perturbative' limit is 
also a 'semi-classical limit' of large k. 

The theory of modular forms should be viewed as a gauge theory 
with a non-abelian and non-compact gauge group- the modular group. 
It has all the essential features of the gauge theories of physics but in 
a much simpler mathematical setting: the group is only countably 
infinite instead of being an infinite dimensional Lie group. Thus there 
is no need for renormalization. By studying modular forms we are 
studying the gauge principle in its purest form without contamination 
by the other complications of quantum field theories. The number 
theory of the last century bears witness to the claim that even this 
simplest of all non-abelian gauge theories is very deep: some of the 
deepest problems of number theory could be solved if we understood 
the spectrum of the Hecke operators. 

In this paper we develop only the analogue of lowest order pertrur- 
bation theory ( 'abelian approximation'). Our other ideas on non- 
abelian gauge theories ('summing planar diagrams') also should have 
analogues here and should lead to deep results in the future. I hope 
an enterprising reader will take up this challenge. 

5.4 Hecke Operators 

We now return to the exposition of the classic theory of modular forms 
due to Hecke. Hecke was motivated by the work of Mordell who in 
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(98) 



turn was trying to understand the Ramanujan conjecture on the r- 
function. 

A lattice on the complex plane is a set of points 

w ~ w + ruj\ + su>2, r,s E Z; (97) 

the fundamental region is a parallelogram with side z. A linear change 
of basis with integer coefficients and determinant onef^] ad — be = 1, 

\UJ2/ \c d) \uj 2 J 
does not change the lattice. By identifying the opposite sides of this 
parallelogram we get a torus. By a rotating our co-ordinate system 
and choosing an appropriate unit of length we can choose U2 = 1. 
Also, we can reflect around this axis if needed to make uj\ lie in the 
upper half plane. Thus only the ratio z = ^ is needed to specify a 
lattice. A modular transformation 

z ^az + b ad-bc=l (99) 
cz + d 

is just the effect of a change of basis on this ratio. We call this lattice 
L z . 

A fractional linear transformation 

az + b 



(100) 



cz + d 

with integer coefficients will map L z to a sublattice of indexF] n = 
ad — be. Each sublattice of index n corresponds to an orbit of T(l) on 



19 If the determinant is not one we change the area of the fundamental domain. See below. 
20 This means that there are n fundamental domains of L z in one fundamental domain 

of the sublattice. 
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the set 

(fa b\ ) 

T n = < I J \a, b, c, d G Z, ad — be = n > (101) 

since an action by T(l) would not have changed the original lattice 
L z . ( For n / 1, T n is not a group; r m r„ C T mn ). 

We now define the action of the Hecke operator T{n) on a modular 
form / as a sum over all the sublattices of indexF 1 ] n: 



[T(n)f]{z) = n*- 1 £ f(h(z)) 
h&r n /r 



dh(z) 



dz 



k/2 



(102) 



Using the fact an action by T(l) merely permutes the terms of this 

sum ( the left action of T(l) on the coset T n /T(l)) we can show that 

T(n)f is also a modular form of weight k. 

By a right action of T(l) we can bring any element of T n to the 

upper triangular form ; actually we can enumerate the elements of the 
fa b\ 

coset r n /r(l) by with ad = n, b = 0, 1, • • • d — 1 . For proofs 



d, 



sec 



18 1 . A more explicit formula for the Hecke operator is thus, 

[nn)/](,) = ^E^E/(^)- (103) 

ad=n b=0 

In particular, for prime p, 

[T(p)f](z)=p k ~ 1 f(pz) + -J2f (—) ■ (104) 



In terms of Fourier coefficients: 

oo oo 

f(z)=Y,f m e 2nimz ,[T(n)f](z)=Y / 7n(m)e 2 ™ mz , (105) 



a b\ 



21 We denote h(z) = ff±^ for h = 



c d I 
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where, 

7n(m)= V d^ 1 /^- (106) 

d| (n,m) 

In particular^], 

[T( P )f] m = fpm + 5{ P \m)p k ~ l f^ (107) 

It follows then that T(mn) = T(m)T(n) if m and n are coprime. 
More generally we can show 

T(m)T(n)= £ d^T (™V (108) 

d|(m,n) 

In particular, the Hecke operators commute with each other. It is 
also useful that for prime powers we have a recursion relation, 

T(f +l ) = T(p)T(p r ) - p^Tif- 1 ) (109) 

which can be solved in terms of Tchebycheff polynomials^: 

T( P r ) = p^u r QpVr(p)) . (in) 

Thus the T(p) for prime p determine all the T(n). 

It is not difficult to establish an inner product on Sk with respect 
to which T(n) are hermitean: 



<h,h>= tt(z)f 2 (z)[Imz] k - 2 d 2 z, (112) 



Id 

22 5(p\m) = 1 if p divides m and zero otherwise. We use Y^,b=o e 5 = p3(p\ m ) to 

derive this formula. 

23 The Tchebycheff polynomials are defined by 

U (x) = 1, Ui(x) = 2x, U r+1 (x) = 2xU r {x) - t/ r _i(x). (110) 
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where, 

D = {z = X + iy\-^<X<^, X 2 +y 2 >l}. (H3) 

The point here is that y k ff(z)f2{z) is a modular invariant function, 
so that it can be integrated after multiplying by the modular invariant 
volume form dx ^f y to get a modular invariant quantity. Of course since 
each tile contributes the same amount we must restrict the integral to 
one fundamental domain. 

Thus we have a set of commuting hermitean matrices, there is an 
orthogonal basis of simultaneous eigenvectors, with real eigenvalues: 
the Hecke forms. 

Suppose we have as simultaneous eigenvector a cuspform satisfy- 
ing, 

T(n)f(z) = X n f(z), Vn. (114) 

The convention is to normalize an eigenvector by setting the first 
Fourier coefficient /i = 1. Taking Fourier coefficients of both sides, 
we get f n = X n : the Fourier coefficients of a simultaneous eigenvec- 
tor of the Hecke operators are the eigenvalues. It follows that these 
coefficients satisfy a multiplicative identity 

fmfn= E d^f^- ( 115 ) 

d\(n,m) 

There is a Dirichlet series associated to any modular form 

OO OO p 

/(z) = / + E/ne 2 ™ 2 , 0(s) = E Z 7 ( 116 ) 

1 1 n 
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or, 



r( s ) y 



y s -V(iy)-Mdy. (ii7) 



The modularity of / implies a functional equation^] for <fi: 

(2tt)~ s i»0(s) = (-l)f (2vr) s - fc r(A; - s)<j>(k - s) (119) 

and conversely. The multiplicative property of the coefficients of the 
Hecke eigenforms yields a product formula: 

^( s )=rii — f - s \ k-i =g - ( i2 °) 
P 1 ~ fpp +p p 

This is reminiscent of the Riemann zeta function C( s )- Indeed the 
theta function, of which £(s) is the Mellin transform, is a modular 
form of a congruence subgroup of level 2. 

It is of much interest to understand the behavior of the eigenvalues 
of the Hecke operators. They have been related to the zeros of the zeta 
function of algebraic varieties over finite fields (Eichler, Sato, Deligne). 
There are many deep conjectures about the behavior of the eigenvalues 
A(n) for large n. The simplest case is when the dimension of the space 
of cusp forms is one: when k = 12 the only cusp form is the function 
A(z) we introduced earlier. In this case the Hecke operators are lxl 
matrices: just numbers. From the above it is clear that these numbers 



24 To see this, just take the Mellin transform of the condition for invariance under 
inversion: 

/(-V = /(*)• (118) 

z 
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are just the Fourier coefficients of the function A(z). In other words, 
for the case k = 12, the Hecke operators reduce to the Ramanujan 
r-function: T(n) = r(n). In fact Hecke discovered these operators by 
generalizing some ideas of Mordell on the modular form A(z) to the 
case of higher weight. 

From our earlier discussion, we are led to consider the limit of large 
level where the invariance group becomes abelian. We now present a 
simple analogue of the Hecke problem for this case of periodic func- 
tions. 

5.5 Hecke Operators on Periodic Functions 

Any modular form is periodic so we can expand it in a Fourier series 
f( z ) = En=i fn&n{z)- Thus it is tempting to think of the space of 
modular forms as a subspace of the space of periodic functions V. To 
make this idea precise, we would like to have a norm on V. 

The inner product on modular forms given above can be written 

as 

/I „ oo 

* dx / dyy k - 2 f(z)f(z) = ]T f*f m 9nm (121) 

where 

(y < 7(1 " * 2 )) dy (122) 

is a positive sesquilinear form. We can extend the integral to the 
fundamental region of the translation group to get an inner product 
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on V: 

/ f OO l ^~ > 
2 dx / y k ~ 2 rHy = ^3mK m . (123) 

In terms of the basis e m , we have h nm as a diagonal sequilinear form 
Km = (en, e m ) = (k — 2)\[Aitn) 1 - k 5 nym . (124) 

The Hecke operators are given by quite simple formulae in terms 
of the Fourier components. We can work out easily their dual action 
on the basis e m (z): 

T(p)e m (z) = p k - l e mp {z) + 5(p\m)e f (z) (125) 

for prime p. This can then be used extend to them as operators on 
V. The spectral problem for T(p) in the infinite dimensional space V 
is much simpler than its counterpart in Sk- We will solve this simpler 
problem and see that it has close connections to the theory of random 
matrices. 

How will we recover the Hecke operators on modular forms? We 
could try thinking of modular forms as a subspace of the space of 
periodic functions. But under the above inner product on V they will 
not be square integrable. The reason is precisely modular invariance: 
each fundamental region contributes an equal amount to the integral, 
and the region — ^ < x < \ contains an infinite number of such 
regions. Thus, the expansion of a modular form in the basis e m (z) is 
not convergent in the above norm (., ). On the other hand, the integral 
for < ., . > corresponding to the sesqilinear form g nrn also can be 
extended to an inner product on V. It is convergent on modular forms 
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but is a degenerate sesquilinear form in V: the integral is restricted 
to the region x 2 + y 2 > 1. We can quotient V by the null space of g nm 
to get a finite dimensional space Sk that is a 'gauge fixed' version of 
Sfc. That is, instead of thinking of modular forms as a subspace of V, 
we think of them as a quotient of V by the null space of g nm . The 
'gauge fixing' amounts to the choice of one particular fundamental 
region among the infinite number as the domain of the integral. Thus 
the spectral problem of the Hecke operators on modular forms can be 
replaced by that on the space Sk- 

As k grows the dimension of Sk grows linearly with k. We will see 
that in a certain sense the two inner products g nm and h nm approach 
each other. Thus the discrete eigenvalues of T(jp) are so close together 
as k — ► oo that they merge to the form the continuous spectrum of 
T(p) on V. 

5.6 Toeplitz Operators 

We now solve the spectral problem for the Hecke operators on periodic 
functions. It will be more transparent to transform to the orthonormal 



There is then a simple description in terms of Toeplitz operators. 
The operators 



basis |m >= (4-irm) 1 k (k — 2)1 2 e m we have 




(126) 



A^(p)\m >= \pm >, A(p)\m >= S(p\m)\ — > 



(127) 



49 



are adjoints of each other and satisfy 

A{p)A\p) = l,A(p)A(p') = A(p)A(p'),A(p)AHp') = A\p')A(pXl28) 

for p ^ p' . That is, they form a commuting set of Toeplitz operators 
labelled by the prime numbers. Being isometries \A\ = \A'\ = 1. It 
follows easily that T{p) are hermitean and that 

\T{p)\<2p h ^r . (129) 

The analogue of this inequality on the space of modular invariant 
functions (rather than periodic functions) is a much deeper statement. 

5.7 Connection to Random Matrices 

The simultaneous eigenfunctions of the T{p) can now be obtained in 
terms of Tchebychev polynomials. The spectrum is connected with 
the Wigner distribution for random matrices. 

It is enough to study each T(p) separately: the theory 'localizes' 
completely. To see this, represent each number m in terms of its prime 
decomposition: 

m = Y[p Up . (130) 

p 

The product is over the set of all primes; v p = 0, 1, ■ ■ ■ with only a 
finite number of them being non-zero. Then 

A\p)\u 2 ,u 3 ,--- >= \u 2 ,--- u p + !,-■■ >, (131) 
A(p)\» 2 , v 3 , • • • >= 5(u p + 0)1^, ■ ■ ■ u p - 1, • • ■ > . (132) 
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Thus A(p),A^(p) act only on the p-th. entry. 

The Toeplitz algebra is the associative algebra generated by a pair 
of elements satisfying the relation 

AA^ = 1. (133) 

The standard representation in terms of an orthonormal basis \v > 
, v = 0, 1 ■ • ■ is 

A*\v>=\v + 1>, A\v >= 5{v ^ 0)\v - I > . (134) 

This is precisely the representation that we have. 

Voiculescu |2(J has found a remarkable connection between the 
theory of random matrices and the Toeplitz algebra. Given any poly- 
nomial f : R —> R, define 

J trf(X)e-^^ x2 dX 
< / >^= J-^x^X (135) 

the integral being over over all hermitean N x N matrices. Thus X is 
a hermitean matrix whose matrix elements are independent random 
variables. Then, Voiculescu shows that 

lim < / >tv=< 0|/ (A + A^) |0 > . (136) 

There is a probability distribution on R, the Wigner semi-circle dis- 
tribution, such that 

lim < / > N = f f(x)p(x)dx. (137) 

Explicitly, 

p(x) = 6(\x\<2)^[A-x 2 ]. (138) 
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Thus we see that the Hecke operators T{p) on periodic functions 

fc-l k-2 

are hermitean operators whose spectrum is the interval [—2p 2 y 2p 2 ]. 
The (generalized) eigenfunctions are given by Tchebycheff polynomi- 
als. The Wigner distribution gives the spectral density. Thus each 
T(p) behaves like a hermitean random matrix; the different Hecke op- 
erators for different prime p commute with each other, so they are not 
'free' in the sense of Voiculescu; instead they are statistically indepen- 
dent in the more conventional sense. 



5.8 The limit of Large Weight 

Recall that the main difference between the exactly solvable model 
above and the the theory of modular forms is that we replaced the 
sesquilinear form g nm by the simpler one h nm . We now show that in 
the limit of large weight k this is a small correction so that what we 
obtained above is the asymptotic behavior as k — > 00. 
Note that we can split 

hnm = 9nm ~\~ Qnrn (1^9) 

where k is the integral over the complimentary region 

2 dxe 2ni{m - n)x / y k - 2 e- 27r(m+n)y dy. (140) 

-i Jo 

If we rewrite this in the orthonormal basis of h nm , we will get 

9nm + Qrim (141) 

where q nm = q nrn / ' y/ (h nn h mm ) etc . We will show that \q nm \ tends to 
zero as k — > oo. 
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Now, evaluating the y-integral, 
l-fc 

9 

Qnm 



dxe 2ni(m~n)x _ 



1 



(fc-2)! 

where the incomplete Gamma function is defined by 



r (k - 1, 2vr(m + n)y/(l - a$42) 



(143) 



We will study this limit as A; — > oo keeping m, n fixedp^|. Now recall 
that as s — 1 > u, the maximum value of the integrand is attained at 
its upper limit in this case, so that T(s,u) < u s e~ u . Then 



< 



rra+ra 
2 

y/(mn) 



l-k 



2ir{m + n)V(l - ^ 2 



(fc-2)! 
fe-i 



Again, replacing the integrand by its largest value ( which is attained 
at x = 0), we get 

[47iy(mn)]* -1 



< 



(145) 



(fc-2)! 

The growth of the factorial beats the exponential growth for fixed m 
and n. 

Thus we can see why the distribution of eigenvalues of Hecke op- 
erators resemble those of random matrices by our semi-classical ap- 
proximation method. Moreover we see why the spectrum of T(p) is in 

fc-i fe— i 
the interval [— 2p 2 } 2p 2 j. 



25 We should really be estimating the operator norm of q. I hope that the arguments 
here will motivate a more rigorous analysis. 
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5.9 The Poisson Algebra of Modular invariant 
functions 

We should expect that in the limit of large weight, the theory of mod- 
ular forms is well-approximated by a classical theory. More precisely 
the algebra of matrices on Sk should tend to the Poisson algebra of 
functions on U/T(l). In particular there will be functions on the upper 
half plane which are classical approximations to the Hecke operators. 
These Hecke functions will have vanishing Poisson brackets relative 
to each other. Their range will give a classical approximation to the 
Hecke eigenvalue problem. This is another, (manifestly gauge invari- 
ant) of studying the limit of large weight. 

The upper half plane is a symplectic manifold with symplectic 
form 



This means that x and \ are canonical conjugates: 



The set of functions on the upper halfplane form a Poisson algebra 
with the bracket 



{y-\x} = \. 



(147) 




(148) 



The modular transformations 



xi— ¥ x + l,yi— ► y 



(149) 
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and 



x 



y 



(150) 



x 



.2 + y. 




x- 



X 



are canonical transformations. Thus the space of modular invariant 
functions is a sub-Poisson algebra. The quotient of this by its center is 
the algebra of 'gauge invariant observables' if we regard the modular 
group as a 'gauge group'. We can construct such observables from 
smooth functions of the upper half plane (vanishing sufficiently fast 
at infinity) by averaging over orbits. The Maas forms provide nice 
examples of such 'observables'. 

In the limit of large k, the Hecke operators should tend to cer- 
tain modular invariant functions ( which are not holomorphic) that 
have zero Poisson brackets relative to each other. The range of these 
functions is the large /c-limit of the Hecke spectrum. We should also 
be able to derive a systematic semi-classical expansion in powers of 
P But this paper is getting long already; I hope to return to these 
questions in a later publication. 
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